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Abstract 

Let us consider the simplest model of one-dimensional probabilistic cellular automata 
(PCA). The cells are indexed by the integers, the alphabet is {0, 1}, and all the cells evolve 
synchronously. The new content of a cell is randomly chosen, independently of the others, 
according to a distribution depending only on the content of the cell itself and of its right 
neighbor. There are necessary and sufficient conditions on the four parameters of such a PCA 
to have a Bernoulli product invariant measure. We study the properties of the random field 
given by the space-time diagram obtained when iterating the PCA starting from its Bernoulli 
product invariant measure. It is a non-trivial random field with very weak dependences and 
nice combinatorial properties. In particular, not only the horizontal lines but also the lines 
in any other direction consist in i.i.d. random variables. We study extensions of the results 
to Markovian invariant measures, and to PCA with larger alphabets and neighborhoods. 
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1 Introduction 

Consider a bi-infmite set of cells indexed by the integers Z, each cell containing a letter from a 
finite alphabet A. The updating is local (each cell updates according to a finite neighborhood) , 
time-synchronous, and space-homogeneous. When the updating is deterministic, we obtain a 
Cellular Automaton (CA), and when it is random, we obtain a Probabilistic Cellular Automaton 
(PCA). Alternatively, a PCA may be viewed as the discrete-time and synchronous counterpart 
of a (finite range) interacting particle system. We refer to [12] for a comprehensive survey of 
the theory of PCA. 

There are two complementary viewpoints on PCA. First, it defines a mapping from the set of 
probability measures on A z into itself. Second, it defines a discrete-time Markov chain on the 
state space A 1 '. A realization of the Markov chain defines a random field on „4 ZxN , called a 
space-time diagram. An invariant measure for a PCA is a probability measure on A 1, which 
is left invariant by the dynamic. Starting from an invariant measure, we obtain a space-time 
diagram which is time-stationary. Our goal is to study the stationary random fields associated 
to some particular and remarkable PCA. 

First, we consider the image by a PCA of a Bernoulli product measure. The resulting measure 
is described via explicit formulas for its finite-dimensional marginals. Second, we use this 
description to revisit a result from pQ (see also |13|. [T2] ) with a new and simple proof: explicit 
conditions on a PCA ensuring that a Bernoulli product measure is invariant. Third, we focus 
on the equilibrium behavior of PCA having such a Bernoulli product invariant measure. The 
resulting space-time diagram turns out to have an original and subtle correlation structure: it is 
non-i.i.d. but, in any direction, the "lines" are i.i.d. In the case of an alphabet of size two and a 
neighborhood of size two (the updating of a cell depends only on itself and its right-neighbor) , 
the stationary space-time diagram satisfies additional remarkable properties: it can also be seen 
as being obtained by iterating a transversal PCA in another direction. 

The paper is structured as follows. General definitions are given in Section 2. A special emphasis 
is put on the simplest non-trivial PCA, that is, the ones defined on an alphabet of size 2 and 
a neighborhood of size 2. They are studied in details in Section 3 and 4. In Section 5, we 
consider the extension to general alphabets and neighborhoods, and we also consider the case 
of Markovian invariant measures. In Section 6, we revisit classical results on CA in view of the 
PCA results. 

Notations. Given a finite set A, the free semigroup generated by A is denoted by A + . The 
length, that is, number of letters, of a word u £ A + is denoted by \u\. The number of occurences 
of the letter a £ A in a word u G A + is denoted by |it| a . 
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2 Probabilistic cellular automata (PC A) 



2.1 Definition of PCA 

Let A be a finite set, called the alphabet, and let X = A z . The set Z will be referred to as 
the set of cells, whereas X is the set of configurations. For some finite subset if of Z, consider 
U = (yk)keK G .A^- The cylinder defined by y is the set 

[y] = {x G X | V/c G if, x fc = yfc} • 

For a given finite subset if, we denote by C(if) the set of all cylinders of base if. Given 
if, L C Z, we define if + L = {u + f|uGif, vel). 

We denote by 7W(»4) the set of probability measures on A. Let us equip X with the product 
topology, which can be described as the topology generated by cylinders. We denote by M.(X) 
the set of probability measures on X for the Borelian cr-algebra. 

Definition 2.1. Given a finite set N C Z, a transition function of neighborhood N is a 
function f : A^ — >■ At (y4). TTie probabilistic cellular automaton (PCA) of transition function f 
is the application F : Ai(X) — > Ai(X), n i — > fiF, defined on cylinders by: Vif, Vy = (yk)k£K, 

[x}eC(K+jKf) k&K 

Assume that the initial measure is concentrated on some configuration x £ X. Then by appli- 
cation of F, the content of the k-th cell is updated to a £ A with probability f((xk+ v )vej\f){0')- 

We keep the notation / for the extended mapping M{A N ) ^M(A), v^vf with 

Va G A, uf(a)= v{u)f{u){a). 

2.2 Space-time diagrams 

A PCA is a Markov chain on the state space X. Consider a realization (X n ) ng N of that Markov 
chain. If X° is distributed according to fj, on X, then X n is distributed according to fiF n . The 
random field (X n ) n( zfq = (X'^)k£Z,neN is called a space-time diagram (the space-coordinate is k, 
and the time-coordinate is n). 

If the neighborhood is M = {0, 1}, for symmetry reasons, a natural choice is to represent the 
space-time diagram on a regular triangular lattice, as in Figure [TJ 

The dependence cone T>(i, n) of the variable Xf is defined as the set of variables which are 
influenced by the value of Xf. If the neig hborhood is M = {0, . . . £}, then V(i, n) = {A"+ fc , k G 
N, -kl < j < 0}. 

Next lemma follows directly from the definition of a PCA. 

Lemma 2.2. Let (i,n) belong to Z x (N \ {0}) and let S be a subset of Z x N such that 
T>(i,n) n S = 0. Then, Xf is independent of {XJ 1 )^^^ conditionally to (X^J) v( zj^. 

We point out that if a PCA has positive rates, i.e., Vn G A^ , Va G A, f(u)(a) > 0, then any of 
its stationary space-time diagram is a Markovian random field. We refer to [10] for an in-depth 
study of the connections between Gibbs states and stationary space-time diagrams of PCA. 
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Figure 1: Space-time diagram. 



2.3 Product form invariant measures 

Definition 2.3. For p € [0,1], we denote by fx p the Bernoulli product measure of parameter p 
on {0, 1} Z , that is, [i p = B® 1 ' , where B p denotes the Bernoulli measure of parameter p on {0, 1}. 
Thus, for any cylinder [x], we have 

Hp[x\ = (l-p)M°pNi • 

We give a first property of the space-time diagram that is shared by every PCA having a 
Bernoulli product invariant measure. 

Lemma 2.4. Let F be a PCA of neighborhood {0, ...,£}. Assume that \x v F = fi p and consider 
the stationary space-time diagram obtained for that invariant measure. Then for any a > 
the line L Q = {(k,n) G Z x N | n = ah} is such that the random variables (X% ){k, n )eL a are 
i.i.d. 

Proof. Let us show that any finite sequence of consecutive random variables on such a line is 
i.i.d. We can assume without loss of generality that the first of these points is Xq. Then, using 
the hypothesis on the slope, we obtain that the other random variables on that line are all 
outside the dependence cone of Xq. Thus, the (n — l)-tuple they constitute is independent of 
Xq. By induction, we get the result. □ 

3 PCA of alphabet and neighborhood of size 2 

For the time being, we assume that the neighborhood is N = {0, 1} and that the alphabet is 
A = {0, 1}. For convenience, we introduce the notations: for x, y 6 A, 

Oxy = 0l y = /(*,!/)(!), 0% = f(x,y)(0) = 1 - 6 xy . 
Observe that a PCA is completely characterized by the four parameters: 8qq, 9q\, 9\q, and On. 

3.1 Computation of the image of a product measure by a PCA 

The goal of this section is to give an explicit description of the measure fJ, p F, where \x v is the 
Bernoulli product measure of parameter p, as a function of the parameters 9qq, #oi> &io, @n- 
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Let us start with an observation. Consider (Y n ) nf zz ~ [J-pF. Clearly, we have: (Y"2n)neZ ~ fJ> p 
and (^2n+i)nez ~ /V But the two i.i.d. sequences have a complex joint correlation structure. 
It makes it non-elementary to describe the finite-dimensional marginals of n p F. 

Assume that the parameters satisfy: 

(0oo,0oi), (0io,0u) 2 {(0,0), (1,1)}. (1) 
For p G (0, 1), a G {0, 1}, define the function 



g a : [0,1] — >• (0,1) 

q ' ^ (1-g) (l-p)0 o t o + (l-9)P^oi + 9(l-p)^o + 9P0 



(2) 



ii 



Consider three random variables Xo, Xi, Y"o with (Xo, Xi) ~ B q ® B p and Yb ~ (B q (g) -Bp)/- In 
words, 50(5) is the probability to have Yq = a. With the condition ([l]), we have g a (q) € (0, 1) 
for all q. Observe also that: go(q) + gi(q) = 1 for all q. 
For p G (0,1), a G {0, 1}, we also define the function 

h a : [0,1] — ► [0,1] 

5 •— > [(l-g)p0& + ?P01lW<O~ 1 - U 

Consider Xq,Xi,Yq with (Xo,Xi) ~ B q ®B p and Yb ~ (£><j &> •Sp)/- In words, /i a (g) is the 
probability to have X\ = 1 conditionally to Yq = a. 

Proposition 3.1. Consider a PCA satisfying Consider p G (0, 1). For ocq ■ ■ ■ a n -\ G A n , 
the probability of the cylinder [a>o ■ ■ ■ a n -i] under n p F is given by: 

n-l 

fj, p F[a ■ ■ • a n _i] =g ao (p) JJ 5^(^-1 (h ai _ 2 (. .. K {p) .. .))) . 

i=l 

By reversing the space-direction, we get an analog proposition for a PCA satisfying the sym- 
metrized condition: (0 OO , 0io), (6*01, 0n) {(0, 0), (1, 1)}. 

Proof. Let us compute recursively the value n p F[ao ■ ■ ■ a n _i]. We set X = X° and Y = X 1 . 
Assuming that X ~ fj, p , by definition, 

H p F[a ] = F(Y = a ) = 9a (p)- 

We can decompose the probability fi p F[aoai] into 

H p F[a Qai ] = F(Y = a , Y x = «i) = P(*i = a x \ Y = a ) P(Y = a ). 

By definition, the conditional law of X\ assuming that Yq = ao is given by Bh ao ( p y So the law 
of (Xi, X2) is Bh a (p) <8> B p and we obtain 

H p F[olqol\\ = g ai (h ao (p)) g ao (p). 

More generally, we have: 

P(Y = a ...Y k = a k ) = F{Y k = a k \Y = a ... Y k . x = a k ^) F{Y = a . . . Yfc-i = a fc _i). 

By induction, the law of X k knowing that Y = a . . . Y^_i = a k-l is s /i afc _ 1 (/»« fc _ 2 (...Ao (p)-))- 
The result follows. □ 
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3.2 Conditions for a product measure to be invariant 



For x G X, denote by 5 X the Dirac probability measure concentrated on the configuration x. 
The probability measure [i\ = 8y/. is invariant for the PCA F if and only if On = 1. Similarly, 
Mo = ^o z ^ s invariant for F if and only if #00 = 0. 

Using Proposition 3.1, we get a necessary and sufficient condition for p p , p G (0, 1), to be an 
invariant measure of F. The result is stated in Theorem 3.2 It already appeared in pQ and 
|12j . but our proof is new and simpler. 



Theorem 3.2. The measure \l v , p G (0, 1), is an invariant measure of the PCA F of parameters 
$00> $01) #10; 011 */ an d only if one of the two following conditions is satisfied: 

(t) (i-p) e o+pe 01 = (i-p) 9 W + P e 11 =p 
(ii) (1 - p) 9oo + P 0ro = (1 - p) 0oi + p Ou = p ■ 

In particular, a PCA has a (non-trivial) Bernoulli product invariant measure if and only if its 
parameters satisfy: 

0oo (1 - On) = 0io (1 - M or 0oo (1 - 0u) = 0oi (1 - 0io) • (4) 

Proof. Let us assume that F satisfies condition (i) for some p G (0, 1). Then, the function g\ 
is given by gi{q) = (1 — q) p + q p = p, and go(q) = I — gi(q) = 1 — p. By Proposition 3.1 we 
have, 

Vq = a •••«„-! G A n , p p F[a] = (1 - pjMopM* = /^[a] . 
So is an invariant measure. 

Now, assume that the PCA F satisfies condition (ii). Let us reverse the space direction, that is, 
let us read the configurations from right to left. The same dynamic is now described by a new 
PCA F defined by the parameters 0oo = 0oo>0oi = 0io>0io = 0oi)0n = 011- So, the new PCA 
satisfies condition {%). According to the above, we have n p F = fi p . Let us reverse the space 
direction, once again. Since the Bernoulli product measure is unchanged, we obtain fj, p F = fj, p . 

Conversely, assume that (i p F = fi p . It follows from Proposition 3.1 that for any value of the a^, 
we must have gi(h an l (h an _ 2 (. . . h ao (p) . . .))) = p. Since g\ is an affine function, there are only 
two possibilities: either g\ is the constant function equal to p; or h an _ x (h an _ 2 (. . . h ao (p) . . .))) = 
p for all values of ozq, . . . , a n -\ G A. 

In the first case, observe that 

9i(q) =q[-(l-p) 0oo - P 0oi + (1 - P) 0io + P n ] + (1 - p) 0oo + P i ■ 
To get: Vq G [0, 1], gi(q) = p, we must have condition (i). 

In the second case, we must have ho(p) = h±(p) = p and g\ (p) = p. Using go(p) = 1 — p and 
91 (p) =P, w e get: 

h (p) = [(l-p)p(l-e 01 )+pp(l-9 11 )}(l- P )- 1 
hi(p) = [(1 -p) P 0oi +PP 0n]p~ 1 = (1 -p) 0oi +P 0n ■ 

The equality hi(p) = p provides the condition (1 — p) 6q\ + p 9n = P- Let us switch to the 
equality /io(p) = P- We have: 

h (p)=p ^ (1-jO(1-0oi)+p(1-0h) = 1-P 
^=> (1 -p) O i +p On = P ■ 

So, we obtain condition (ii). □ 
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To complete Theorem 3.2 let us quote a result from |13j . We recall that a PCA has positive 



rates if: Vu <E W, Va € A, f(u)(a) > 0. 

Proposition 3.3. Consider a positive-rates PCA F satisfying condition (i) or (ii), for some 
p G (0, 1). Then F is ergodic, that is, [i p is the unique invariant measure of F and for all initial 
measure fi, the sequence (fiF n ) n >o converges weakly to fj, p . 

Assessing the ergodicity of a PCA is a difficult problem, which is algorithmically undecidable in 
general, see [lit [3] . On the other hand, a long standing conjecture had been that any PCA with 
positive rates is ergodic. However, in 2001, Gacs disproved the conjecture by exhibiting a very 
complex counter-example with several invariant measures [5j. In this complicated landscape, 
Proposition |3.3| gives a restricted setting in which ergodicity can be proven. 



Observe that Proposition 3.3 is not true without the positive-rates assumption. Consider for 
instance the PCA defined by: 9 Q0 = p/(l - p), 6 m = 0, 9 W = 0, 6 U = 1 for some p £ (0, 1/2]. It 
satisfies (i) and (ii), but it is not ergodic since 8]% and /i p are both invariant. 



3.3 Transversal PCA 

We assume that fi p is invariant under the action of the PCA, and we focus on the correlation 
structure of the space-time diagram obtained when the initial measure is fi p . Observe that 
this space-time diagram is both space-stationary and time-stationary. By time-stationarity, 
the space-time diagram can be extended from Z x N to Z 2 . From now on, we work with this 
extension. 

Let (X k , 

n)k,neZxZ be a realization of the stationary space-time diagram. 




It is convenient to define the three vectors u, v, and w as in the figure above. The PCA generating 
the space-time diagram is the PCA of direction u. In some cases, the space-time diagram when 
rotated by an angle of 27r/3 (resp. — 2-7r/3) still has the correlation structure of a space-time 
diagram generated by a PCA of neighborhood {0, 1}. In this case, we say that, in the original 
space-time diagram, there is a transversal PCA of direction v (resp. w). 

Proposition 3.4. Under condition (i), each line of angle tt/3 of the space-time diagram is 
distributed according to \x p . Moreover, their correlations are the ones of a transversal PCA of 
direction v and rates given by: $00 = #00 > $01 = #io> $10 = $01; $11 = #11 ■ 



To prove Prop. 3.4, we need two preliminary lemmas. Set X 
have in particular (X, Y) ~ fJ, p F). 



X° and Y = X 1 , so that we 
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Lemma 3.5. Under condition (i), the variables (Yk)k>o are independent of Xq, that is, for any 
n>0, 

n 

P(X = x , (Yi) <i< n = (yi)o<i< n ) = H P [xo] H n P [yi] . 

i=0 

Proof. The left-hand side can be decomposed into: 

^ P((Aj)o<i<n+l = (xi)o<i< n +l, (Yi)o<i<n = (j/i)o<i<n), 

x 1 -x n+1 e{o,i}"+ 1 

which can be expressed with the transition rates of the PCA as follows: 

n 

a;i-2; Il+ ie{0,l}"+ 1 «=0 



H6{0,1} x 2 e{0,l} a;„ + ie{0,l} 



Condition (z) can be rewritten as: 

Va,6,cG {0,1}, Y VpWab = V P [c\- 
fee{o,i} 

Using this, and simplifying from the right to the left, we obtain: // p [xo] 111=0 /-VM- ^ 
Lemma 3.6. Under condition (i), for any n > 0, 

n 

¥(X = X ,Xt = Xl, (Y;) <j<n = (yi)o<i<n) = Mp^ojAipM^lU II VplVi) • 

i=l 

Proof. The proof is analogous. We decompose the left-hand side into: 

^ ¥((Xi)o<i<n+l = (xi)o<i< n +l, {Yi)o<i<n = (Ui)o<i<n), 

x 2 -x n+1 G{0,l} n 

which can be expressed with the transition rates of the PCA as follows: 

n 

y ^p[ x °] n M^+iMk+i 

:E2-z n+ ie{0,l} n i=0 

= n P [xo\n P [xi}ef oxl Y /%>M^U . . . Y ^pi X n+l\ G xZx n+1 - 
x 2 e{o,i} x n+ ie{o,i} 

Using (i) and simplifying from the right to the left, we get the result. □ 



Proof of Proposition \3.4\ To prove the first part of the proposition, it is sufficient to prove 
that the sequence [Xq)^i is i.i.d. For a given n 6 N and a sequence (afc)o<fc<n, let us prove 
recursively that P((Xq ) <fc< n = {ctk)o<k<n) = [ a o • ■ ■ oi n ] . For n = 0, the result is straight- 
forward; and for n = 1, it is a direct consequence of Lemma 3.5 For larger values of n, set 
A = P((X^) < fc < n = (a fc ) <fc<n), we have: 

A = Yj P((^o)o<fc<n. = (Q!fe)o<Jfc<TU (Yi)l<i<n-1 = (Ui)l<i<n-l) ■ 

yi-yn-ieio,!}™- 1 
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Since X$ = Xq, Xq = Yq, it can be rewritten as: 
A = 

yi 

X P(.Yn ■ 

rk 

L o 



= Xq, Xq = Yq, it can be rewritten as: 

F((X^) 2 <k<n = (a k h<k<n | X = a ,Y = Qi, (Yi) 

-yn-i&iO,!}"- 1 



l<i<n-l — {Ui)l<i<n-l) 
i)x<i<n-l = (lli)l<i<n-l) 



(w) 

get: 



i4 = 



o = ao,Y = a\, (Yi) 

iditionally to (Xq, ( Y)o< i<n-i) is equal to the law of (X t 
o<i<n-i- Also, using Lemma 3^ we have: F(X Q = a ,Y = a\, (Xi)l<i<n-i = 
jUp[«o] P(5o = «i, (^.)i<i<n-i = (y«)i<j<n-i)- Coupling these two points, we 

H( X oh<k<n=(pt k ) 2 <k<n\Yo = a 1 , 



, (Yi)l<i< n -l — (Ui)l<i<n-l) 



x n p [a ] P(Y = ai, (Vi)i<i<n-i = (j/i)i<i<n-i" 



= //p[a ] P((A"o)i<fc<„ = (a k )i< k <n) 
By induction, we obtain the result. 




The second part of the proposition consists in proving that 

n 

¥((X^) <k<n = (Pk)o<k<n I (A^) <fc<„+1 = (a k ) <k<n+l) = ]J " ( 5 ) 

k=0 

We prove the result recursively. For n = 0, set A = ¥(X\ = (3q \ Yq = ati, Xq = ceo). We want 
to prove that A = $a° ao . Using the first part of the proposition, we have: 

A = F(Y = ai \X = a Q ,X 1 = Po) F(Xq = a , X x = ft) F(X = a , Yq = a^ 1 

= Oo ^[ a °] MpIAj] MpM" 1 UpWiY 1 = Oalfo tolPd MpM" 1 • 

If a.\ = A) = u, we get A = 6™ ou = i?^ a . Assume that ot\ ^ /3q. Condition (i) can be rewritten 
as: 

MftlOo+^Kld^hl' (6) 

Dividing by // p [cti], we get: 

A = C b q -MA)] M^r 1 = 1 " Cat = = ^ao • 
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For larger n, it is convenient to prove next equality, which is equivalent to ([5]): 

n 

P((^o)o<fc<n+l = («fc)o<fc<n+li (^l)o</c<n = (Pk)o<k<n) = Hp[ a n+l] Hp [ a k]°&^t k+1 a k ■ 

k=0 

The left-hand side can be decomposed into: 

^ P((Xq )o<fc<n+l = («fc)o<fc<n+l, (-^1 )o<fc<n = (Pk)o<k<n, {Yi)2<i<n = (lli)2<i<n) ■ 

V2-y n e{0A}"- 1 

Let us decompose each term of the sum, conditioning by the values of Xq, X\,Yq, and Y\. We 
have: 

P((-^0 )2<fe<n+l = (Q!*)a<k<n+l! (-Xl )2<fc<n = (/?fc)2<fc<n | (^M3, -Xl, ^0, ^l) = ( Q 0, A), a li A)> (^i)2<i<n = (Vi)2<i<n) 
= P((Xq )2<fc<n+l = (Qfc)2<fc<n+1, (Xl )2<fe<n = {Pk)2<k<n \ (lo,Yi) = («1 , /3\ ) , (Yj )2<i<n = {yi)2<i<n) ■ 



and using Lemma 



3.6 



and the equality Hp\Po]^ p = Hpi^i]^^ (see (|6|)): 



P((X ,X 1 ,FoTi) = (a ,/3o,ai,/3 1 ),(y i ) 2 <i<n = {yi)2^<n) 

= Hp[ao]H P [/3o}6T Po F ( Yl = 1 ' ( y *)2<i<n = (j/i)2<i<n) 
= M«oW[ a l]#ai<* = A, (>i)2<i<n = (Vi)2<i<n) 

= A*p[ao]<oo F (( y 0, *i) = ("1, A), (^) 2 <i<n = ( yi ) 2 <i<n) • 

Assembling the pieces together, we obtain: 

)o<fc<n+i = (afc)o<fc<n+i) (-X'l 

= ^["oj^aiao P((A~q )l<fc< n +l = ( a k)l<k<n+l, (^1 )l<fe<n = (/3fc)l<fc<n) • 

We conclude the proof by induction. □ 

Corollary 3.7. Under condition (i), all the lines of the space-time diagram except possibly 
those of angle 2ir/3 consist of i.i.d. random variables. 



Proof. The previous proposition claims that the lines of angle ir/3 are i.i.d. Lemma 2.4 provides 
the result for the lines of angle in [0, vr/3) U (2n/3, tt]. The angles in (tt/3,2tt/3) correspond to 
lines that are outside the dependence cones of the transversal PCA, so we obtain the result by 



applying again Lemma |2.4| for the transversal PCA. □ 
In the same way, one can prove the following. 

Proposition 3.8. Under condition (ii), the lines of angle 2tt/3 of the space-time diagram are 
distributed according to \x v and their correlations are those of a transversal PCA of direction w 
and rates given by ??oo = #oo, #11 = #n and i9 i = 010, #10 = #oi- 

Corollary 3.9. Under condition (ii), all the lines of the space-time diagram except possibly the 
ones of angle it/3 consist of i.i.d. random variables. 

For a PCA satisfying (i) (resp. (ii)), the lines of angle 27r/3 (resp. n/3) are not i.i.d., except if 
the PCA also satisfies condition (ii) (resp. (if). The distribution of the lines of angle 27r/3 (resp. 
7r/3) does not necessary have a Markovian form either. For example, if #oo = 0oi = 1/2 and 
^10 = 0,011 = 1 (condition (i) is satisfied with p = 1/2), one can check that P(A~q = 0, X}l± = 
0,X^ 2 = 0) = 19/64 which is different F(X$ = 0)P(X£ 1 = | X$ = 0)P(X£ 2 = | Xi x = 0) = 
(l/2)(3/4) 2 . 

It is an open problem to know if under condition (i) (resp. (ii)), it is possible to give an explicit 
description of the distribution of the lines of angle 27r/3 (resp. 7r/3). 
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4 Non-i.i.d. random field with every line i.i.d. 



We now concentrate on PCA satisfying both conditions (i) and (ii) for some p G (0, 1). We con- 
sider the stationary space-time diagram associated with (j, p , and we still denote it by (X^)f. tnE z- 

4.1 All the lines are i.i.d. 

For a given p G (0, 1), conditions (i) and (ii) are both satisfied if and only if: 
" 2p — 1 p 



3s G 



P ' 1 ~P 



_ p{l - s) (l-p)a ry x 
oo — — ; j c'oi — "io — s, f 11 — i • (I) 

1-p p 



Example 4.1. For any value of p G (0, 1), the choice s = p is allowed. In that case, the 
transition rates 9ij are all equal to p and the stationary random field is i.i.d., there is no 
dependence in the space-time diagram. 

Example 4.2. If p = 1/2, every choice of s G [0, 1] is valid and the corresponding PCA has the 
transition function f(x,y) = s S x+yuiod2 + (1 - s) 5 x+y+ i mod2 . 



An example of space-time diagram for p = 1/2 and s = 3/4 

Example 4.3. For any value of p G (0, 1/2], it is possible to set s = and then, 6>oi = 8\o = 
0, #n = 1, and #oo = p/0- ~ P)- This PCA forbids the elementary triangles pointing up that 
have exactly one vertex labeled by a 0. 



An example of space-time diagram for p = 1/3 and s = 



Next proposition is a direct consequence of Corollaries 3.7 and|3.9 



Proposition 4.4. Consider a PCA satisfying |?p. Every line of the stationary space-time 
diagram consists of i.i.d. random variables. In particular, any two different variables are 
independent. 



11 



4.2 Equilateral triangles pointing up are correlated 

We have seen that all the lines of the space-time diagram are i.i.d. But the whole space-time 
diagram is i.i.d. if and only if s = p. Indeed, if s ^ p, the random variable X^ +1 is not 
independent of (X£ , X£ +1 ); in words, the three variables of an elementary triangle pointing 
up are correlated. Precisely, the triple (Xj£, X£ +1 , X'^ +l ) consists of random variables which 
are: (1) identically distributed; (2) pairwise independent; (3) globally dependent if s ^ p. The 
"converse" holds. 

Proposition 4.5. Let u be a law on {0, l} 3 such that the three marginals on {0, l} 2 are i.i.d. 
Assume that v is non degenerated (v ^ ^ooO)^ ^ ^lil)- Then v can be realized as the law of 
an "elementary triangle pointing up" in the stationary space-time diagram of exactly one PC A 
satisfying 

Proof. Consider (Xq,X\,Yq) ~ v. Assume that the common law of Xq,X\, and Yq is B p . By 
the pairwise independence, we have: 

p(x = i,x 1 = o ) y = o) = p(x 1 = o,y = o)-p(x = o,x 1 = o,y = o) 

= (i-pf -f(x = o,x 1 = o,y = o). 

We obtain: 

F(X = l,X l = 0,Y = 0) = F{X = 0,Xi = l, y = o)=P(x = o,Xi = o,y = i) 
p(x = o,x 1 = i,y = l) = P(x = i,Xi = o,y = 1) = p(x = i,x 1 = i,y = o) . 

Set q = F(X = 1, Xx = 0, Y = 0) and q x = P(X = 0, Xt = 1, Y Q = 1). We have: 

p(x = o, x x = o, y = o) = (i-p) 2 -qo, nx = i,x x = i,y = i)= P 2 - qi . 

Furthermore: 

go+gi = p(a = o,x x = o,y = i)+p(x = i,Xi = o,y = i) = f(Xi = o,y = i) = p(i-p). 

Using the above, and expressing everything as a function of p and q±, we get: 

P(y = l|X = 0,X 1 =0) = (p(l-p)-qi)/(l-p) 2 
P(y = l|X = 0,X 1 = l) = qxMl-p)) 

p(y = i|x = i,x 1 = o) = <&/(p(i-j>)) 

P(y = l|X = l,X 1 = l) = l-qi/p 2 . 
By setting 9{j = F(Yq = 1 \ Xq = i, X\ = j) and s = qi/(p{l — p)), we recover exactly ([7]). □ 

Proposition 4.6. Consider a PCA satisfying ^ with s ^ p. The correlations between three 
random variables that form an equilateral triangle pointing up decrease exponentially in function 
of the size of the triangle. 

Proof. Let us consider the random field (X|^)fe iTie z. Observe that all its random variables 
are distributed according to B p , and that each line consists of i.i.d. random variables. More- 
over, for any a < b, the variables {X 2 £ +2 ) a <k<b are independent conditionally to the variables 
(-X^) a <fc<&+1- Thus, this "extracted" random field corresponds to the space-time diagram of a 
new PCA, having a neigborhood of size 2 and satisfying Q for the same value of p. To know 
its transition rates 6^ = P(Aq = 1 | X$ = i,X® = j), it is enough to compute 9^ = 6^ . We 
denote this value by 4>(s), since it is a function of s = 9qi = 9io- 
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Summing on all possible values of X° , Xq , X\ (we first consider the case X® = 1 and then the 
one X® = 0), we get: 

0(a) = P [0oi0ii % + (1 - 0oi)0n 0oi + 0oi (1 - 0n) 0io + (1 - 0oi)(l - 0ii) 0oo] 

+(1 -p) [0 O o#oi 0ii + (1 - 0oo)0oi 0oi + 0oo (1 - 0oi) 0io + (1 - OO )(1 - 0oi ) 0oo]- 

Replacing the coefficients 9ij by their expression in function of p and s and simplifying the 
result, we obtain: 

(s - p) 3 



p + 



p(l-p)' 



We proceed similarly for the random field (X^l2)k,neZ- The coefficient 0^ = F(Xq Z = 1 | Xq = 
0,X^ = 1) is equal to </>*(s), which satisfies: 

(p(i-p))V WpO^-pV 

Similar computations can be performed for equilateral triangles pointing up of other sizes. The 
decay of correlation for equilateral triangles pointing up is exponential in function of their 
size. □ 

Next lemma will allow us to characterize completely the triples of random variables that are 
not independent. 

Lemma 4.7. Consider a PC A satisfying Q). The variable Xq is independent of (X^) k ^ . n eN\{o}- 

Proof. Set X = X° and Y = X . It is sufficient to prove that Xq is independent of (ifc)fcgz- 
But (Yk)k>o and (Yk)k<o are independent conditionally to Xq, so that we can conclude with 
Lemma |3.5| and its analog for condition (ii). □ 



Proposition 4.8. Consider a PCA satisfying with s ^ p. Three random variables of the 
stationary space-time diagram are correlated if and only if they form an equilateral triangle 
pointing up. 

Proof. Three variables that form an equilateral triangle pointing up are correlated, see the proof 



of Proposition 4.6. Let us now consider three variables (Z\, Z2, Z3) that do not constitute such 
a triangle. Then, if we consider the smallest equilateral triangle pointing up that contains them, 
there is an edge of that triangle that contains exactly one of these variables. By rotation of 
angle 2ir/3 or translation of the diagram, one can assume that this edge is the horizontal one 



and that it contains the variable Z\, and not the variables Z2,Z%. Now, using Lemma 4.7 
we obtain that Z\ is independent of (Z2, Z3). But since Z2 and Z3 are independent, the three 
variables (Zi, Z2, Z3) are independent. □ 

There are subsets of four variables that do not contain equilateral triangles pointing up and 
that are correlated. It is the case in general of (Xq, X2, Yq, Y\). Let us consider for instance 



the PCA of Example 4.3 The event (Xq, X2, Yq, Yx) = (0,1,1,1) has probability zero, since 
whatever the value of X\ , the space-time diagram would have an elementary triangle pointing 
up with exactly one zero. 




x 2 = 1 
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4.3 Incremental construction of the random field 



Let us show how to construct incrementally the stationary space-time diagram of a PCA satis- 
fying conditions (i) and (ii), using two elementary operations. 

Consider a PCA satisfying (i) and (ii) for some p G (0, 1). Let S C 1? be the finite set of points 
of the space-time diagram that has been constructed at some step. Initially S = {(0,0)} and 

• If (i,n), (i + l,n) G S, (i,n+ 1) G" S, and V(i,n+ 1) n S = 0. Choose X, n+1 knowing 
(X?,X? +1 ) according to the law of the PCA. 

If (i, n), (i, n + 1) G S, (i + 1, n) $ S, and if no point of the dependence cone of (i + 1, n) 
with respect to the transversal PCA of direction v belongs to S: choose Xf +1 knowing 
(X 4 n+1 , A"™) according to the law of the transversal PCA of direction v. 

If (i,n + 1), (i + l,n) G S, (i,n) S, and if no point of the dependence cone of (i,n) 
with respect to the transversal PCA of direction w belongs to S: choose X" knowing 
(X™ +1 , X™ +1 ) according to the law of the transversal PCA of direction w. 

• If (i,n) G" S, and if (j,m) G S implies m < n: choose Xf according to B p and indepen- 
dently of the variables XJ 1 , (J, in) G S. 

If (i, n) S, and if (j, m) G S implies j < i: choose Xf according to B p and independently 
of the variables X™, (j, m) G S. 

If (i, n) S, and if (j, m) G S implies j + m > i + n: choose Xf according to B p and 
independently of the variables X™, (j, m) G S. 

By applying the above rules in the order illustrated by the figure below, one can progressively 
build the stationary space-time diagram of the PCA. Indeed the rules enlarge S in such a way 
that, at each step, the variables of S have the same distribution as the corresponding finite- 



dimensional marginal of the stationary space-time diagram. This is proved by Lemmas 2.2 and 
1371 

On the figure, the labelling of the nodes corresponds to the step at which the corresponding 
variable is computed (after the three variables of the grey triangle). An arrow pointing to a 
variable means that it has been constructed according to the PCA of the direction of the arrow 
(first rule) . The nodes labelled by II are the ones which have been constructed by independence 
(second rule). 
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5 Extensions 



We consider two types of extensions. First, PC A with an alphabet and neighborhood of size 2 
but having a Markovian invariant measure. Second, PCA having a Bernoulli product invariant 
measure but with a general alphabet and neighborhood. 



5.1 Markovian invariant measures 

Markovian measures are a natural extension of Benoulli product measures. In a nutshell, the 
tools of Section [3] can be extended to find conditions for having a Markovian invariant measure, 
but the spatial properties presented in Section [4] do not remain. 

Definition 5.1. Consider a, b S (0, 1). The Markovian measure on {0, 1} Z of transition matrix 

_ / 1 — a a 
Q ~\l-b b 

is the measure vq defined on cylinders by: 

\/x = X m ---X n , VQ[x]=-K Xm Y\ Qx u x 



n-1 

n +1 ' 



where ir = (iro,ni) is such that irQ = ir, -kq + tt\ = 1, that is, tto = (1 — b)/(l — b + a) and 
7ri = a/(l — b + a). 

The Markovian measure vq is space-stationary. If ci — b, then uq — /i a , the Bernoulli product 
measure of parameter a. 

Let us fix the PCA, that is, the parameters (000) #oi, #10, #11) and assume that ([!]) holds. Let us 
fix the parameters a and b in (0, 1) (defining Q and tt as in Definition 5.1). We introduce the 
analogs of the functions defined in ^ and ^ . 

For a £ {0, 1}, define the function: 

g a : [0,1] — > (0,1) 

r 1— > (1 - r) (1 - a) % + (1 - r) a Ofa + r (1 - b) 9f + rb 0ft . { } 

In words, g a {r) is the probability that Yq = a if the law of (Xq, X\) is given by P(Xo = xo, X\ = 
x i) = r x Qx ,xi with ro = 1 — r and ri = r. With condition ([!]) on the parameters, we have 
g a (r) £ (0, 1) for all r. Observe also that: go(r) + 51 (r) = 1. 

For a € {0, 1}, we also define the function: 

fr a : [0,1] — ► [0,1] 

r I—). [(1-r) ae^+rbOf^gair)- 1 . { ) 

In words, h a (r) is the probability to have X\ = 1 conditionally to lo = a if (Xo,Xi) is 
distributed according to the above law. 

Proposition 5.2. Consider the Markovian measure vq and the PCA F as above. For 
olq • • • a n -i G A n , the probability of the cylinder [a± ■ ■ ■ a n ] under vqF is given by: 

n-1 

v Q F[a ■ ■ ■ an-!} =g ao (ir 1 ) g^h^ {h ai _ 2 (. .. h ao (iri) .. .))) . 

i=i 
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Using Proposition |5.2[ we obtain sufficient conditions for having a Markovian invariant measure. 
This provides a new proof of a result mentioned in [T2] and first published in pQ . 

Theorem 5.3. A PCA has a Markovian invariant measure if its parameters satisfy: 

Ooo On (1 - O i) (1 - do) = 0oi Oio (1 - M (1 - 0n) , (10) 
and (8 00 , 6 i), (0 W , On) {(0, 0), (1, 1)} or (0 OO , W ), (0 i,0 n ) ? {(0, 0), (1, 1)}. 



Proof. We treat the case (f9 00 , Ooi), (Oio, On) £ {(0, 0), (1, 1)} (for which Prop. [5^2] holds). The 
case (000) 0io)> (0oi> 0n) {(0, 0), (1, 1)} can be treated by reversing the space-direction. 
Let us assume that the following conditions are satisfied: 



1. for a G {0, 1}, g a (^i) = n a ; 

2. for a G {0, 1}, there exists c a G [0, 1] such that: Vr, h a (r) 

3. for a, (3 G {0, 1}, gp(c a ) = Q a ^. 



Ca] 



Then, by a direct application of Proposition 5.2 the measure vq is invariant. When are these 
conditions fulfilled? 

For a = 1, condition 2 tells us that there exists c\ G [0, 1] such that for any r G [0, 1], 

(1-r) adoi + rbOn =ci((l-r) (1 - a) 6 00 + (1 - r) a 9 01 + r (1 - b) 6 w + rb9n). 
This is the case if and only if: 

a 0oi = ci((l - a) OO + a O i), 6 U = ci((l - 6) iO + & 0n) • 
Thus, condition 2 for a = 1 is equivalent to: 

a (1 - b) 0oi 0io = (1 - a) b #oo 0n • 
In the same way, condition 2 for a = is equivalent to: 

a (1 - b) (1 - 0oi ) (1 - 0io) = (1 - a) b (1 - OO ) (1 - 0n) ■ 



(11) 



(12) 



Eliminating a and b in ( 11 ) and ( 12 ), we obtain the relation (h0| for the parameters of the PCA. 



Conversely, let us assume that relation (10) holds. We will prove that there exists a, b G (0, 1) 
such that the three above conditions are satisfied. 



First observe that ( 11 ) holds if and only if ( 12 ) holds. So, we have a first relation to be satisfied 

nth: 

(13) 



by the parameters a, b G (0, 1) which is ([11 ). Under this relation, condition 2 is satisfied with: 

a(l-%) 6(l-0n) 



(1 - a) (1 - OO ) + a (1 - 0oi ) (1 - 6) (1 - &w) + b(l- 9 U ) 



and 



c\ 



a O i 



bOn 



(1 - a) OO + a 0oi (1 - b) 0io + b On ' 
Now consider condition 3 for a = /3 = 1. Symplifying using ( |14[ ), we obtain: 

Si(ci) = Qn = b ^ (1- a) oo = 6(l-0ii). 
Condition 3 for other values of a and /? provide the same relation after simplification. 



(14) 



(15) 
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Let us show that if equations (11) and (15) are satisfied, then the PCA also fulfills condition 1. 



Is is sufficient to prove that gi(iri) 
(11), we obtain: 

o(l - b) (1 - O i 
Applying the definition Q, we have: 
1 



7Q. Expanding both sides of (12) and simplifying using 
a)b (1 



'10 



) = (1 



7 00 



(16) 



l-b + a 



(l-6)(l-a) 0oo + (1 



b)a 6»oi + a(l - b) Q w + a& 6>n 
_ b ) _ (i _ a ) b (i _ 0OO _ With (1151), 



Using (16), we can replace a(l — b)(0oi + 0io) by a(l 
we finally obtain 5i(vri) = a/{l — b + a) = tt\. 
Now, observe that the system: 

(1 - a) b <?oi 0io = a (1 - b) 0oo 0n 
(I -a) oo = 6(l-0 n ) 

has a unique solution (a, 6) 6 (0, l) 2 . Let Q be the matrix associated with (a,b). Since the 
three above conditions are satisfied, the Markovian measure vq is invariant by the PCA. □ 



(17) 



In the Markovian case, unlike the Bernoulli case, there is no simple description of the law 
of other lines in the stationary space-time diagram. Nevertheless, the stationary space-time 
diagram has a different but still remakable property: it is time-reversible, meaning it has the 
same distribution if we reverse the direction of time. This is proved in |13j . 

Bernoulli product measures are special cases of Markovian measures. Therefore it is natural 
to ask whether all the cases covered by Theorem 3.2 are retrieved in (10). The answer is no. 
Indeed, the measure vq is a Bernoulli product measure iff a = b. Simplifying in (17) and (10), 
we obtain: 

[0oo = 0oi)0n = 0io] or [0oo = 0iO)0n = 0oi] • 
The corresponding PCA have a neighborhood of size 1. This is far from exhausting the PCA 
with a Bernoulli product measure. 



Finite set of cells. It is also interesting to draw a parallel between the result of Theorem 



5.3 and Proposition 4.6 of Bousquet-Melou [2]. In this last article, the author studies PCA 
of alphabet A = {0,1} and neighborhood M = {0,1}, but defined on a finite ring of size 
N (periodic boundary conditions: Xn = Xq), and proves that the invariant measure has a 
Markovian form if the parameters satisfy the same relation (10) as in the infinite case. The 
expression of the measure is then given by: 



N-l 



F(X = x , Xi = xi,..., Aiv_i = xjv-i) 



In 



where Z is a normalizing constant, and where the coefficients o and b defining the matrix Q are 



the solution of the same system (17) as in the infinite case. 



For a PCA satisfying condition (10), we have a Markovian invariant measure both on a finite 



ring and on Z. This is not the case for Bernoulli product measures: except when the actual 



neighborhood is of size 1, PCA satisfying the conditions of Theorem 3.2 do not have a product 
form invariant measure on finite rings. 

Example 5.4. Consider for instance the PCA of transition function f(x, y) = (3/4) 5 x + y mo d 2 + 
(1/4) 5 x +y+imod2 (Example 4.2), on the ring of size 4. Its invariant measure fi is different from 
the uniform measure: 

^(0000) = 573/8192, /i(0001) = 963/16384, //(0011) = 33/512, 
jt(OlOl) = 69/1024, jt(Olll) = 957/16384, /x(llll) = 563/8192 . 
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5.2 General alphabet and neighborhood 

In this section, the neighborhood is M = {0, . . . ,£} and the alphabet is A = {0, . . . , n}. For 
p = (po, . . . ,p n ) such that po + • ■ ■ + p n = 1, we still denote by \i v the corresponding Bernoulli 
product measure on A z . 

For convenience, we introduce the following notations: Vxo, . . . ,xi G A, V7c G A, 

d x -x e = f(x ,...,x e ){k) . 

We define new functions and h^, which generalize the ones in ^ and (|3j). These new functions 
gk and /ifc are not functions of a single variable, but of probability measures on A . Assume 
that: 

VA; G A, Vx ■ ■ ■ xe-i G A e , 3i G A, Ol -x t _ t i > • ( 18 ) 

Let us define: 

g k :M(A e ) — )• (0,1), 

V i — > the probability that Y = k if (X , . . . , X t ) ~ 2? <g> B p , 



h k :M(A e ) 
V 



M(A e ), 

the distribution of (X\, . . . , Xg) conditionally to Yq = k 
if (X Q ,...,X t ) ~V®B p . 



We have the following analog of Proposition |3.1 



Proposition 5.5. Consider a PCA satisfying (18). Consider p = (pi)i£A with pi > for all i. 
For ao ■ ■ ■ a n -\ G A n , the probability of the cylinder [ao • • • «n-i] under fi p F is given by: 



n-l 



(J,pF[a ■ ••«„-!] = g ao (B® £+1 ) Y[ ffa l (^a l _i(^a l _ 2 (- ■ • h ao (B. 



v > 



•))) 



By reversing the space-direction, we get an analog of Proposition 5.5 under the symmetric 
condition: V7c G A,Vx ■ ■ ■ x e _i G A £ , 3i G A, 0^,, x > 0. 



Applying Proposition 5.5 we otain the following result. It already appears in [13] in a more 
complicated setting. 

Theorem 5.6. Consider p = (pi)i£j\_ with pi > for all i. The measure fi p is an invariant 
measure of the PCA F if one of the two following conditions is satisfied: 



Vxq, • • • ,X£-i G A, V7c G A, 



i v x ~-x e _ii 



n 



Vxq, • • • ,X£-i G A, V7c G A, y^^Pi 



ixq---X£_i 



Pk 



Pk 



(19) 
(20) 



Proof. Let us assume that F satisfies condition (19). Then, the function g^ is constant. Indeed, 
9k{V) = 



Pk 



i&A,X( i ---xi_ 1 &A l 



By Proposition 5.5, we obtain that fi p F = [i p . 

Now, like in the proof of Theorem |3.2[ w e can reverse the space direction and define a new PCA 
F. The PCA F satisfies condition (20) iff the PCA F satisfies condition (19). Therefore, if F 
satisfies condition (20), then we have fi p F = fi p , which implies in turn that \x p F = jx p . □ 
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As opposed to Theorem |3.2[ the conditions in Theorem |5.6| are sufficient but not necessary. To 
illustrate this fact, the simplest examples are provided by PCA that do not depend on all the 
elements of their neighborhood. Consider for instance the PCA of alphabet A = {0, 1} and 
neighborhood M = {0,1,2}, defined, for some a,b G (0,1), by: Vit, v G A.^9\ Qv = a, 9^ lv = b. 
This PCA has a Bernoulli invariant measure, but if a ^ b, it satisfies neither condition (19), 
nor condition (20). 

Let us state a result from [13], which extends Proposition 3.3, and completes Theorem |5.6| (For 



the relevance of this result, see the discussion following Proposition 3.3 



Proposition 5.7. Consider a positive-rates PCA F satisfying condition (19) or (20), for some 
P = (,Pi)ieA> Pi > for all i. Then F is ergodic, that is, fi p is the unique invariant measure of 
F and for all initial measure ji, the sequence (/^i ?ri )n>o converges weakly to \i p . 



Condition (|19j) implies that the variables Xq, . . . , X^_i, Yq are mutually independent, since for 
any v G {0, l} e and a G {0, 1}, we have P((A , . . . , X^ x ) = v, Yq = a) = /jl p [v] J2ieAP* ®vi = 
Hp[v]fi p [a]. Similarly, condition (20) implies that the variables X\, . . . , Xi, Yq are mutually 
independent. 



Next lemma is a generalization of Lemma 4.7 



Lemma 5.8. Under conditions (19) and (20), the variable Xq is independent o/(A^)fc g 2 )n gN\{o}- 



Proof. Set X = X° and Y = X 1 . Like in Lemma 4.7 it is sufficient to prove that Xq is 
independent of Y = (Yk)k&i- Let us fix some a,b G Z, (a < < 6), and prove that Xq is 
independent of (Y a , Y a+ i, . . . , Yj,). We have: 



E 



S = P(X = Xq, {Yi) a <i< h = (yi) a <i<b) 

a<i<b+£ = (Xi)a<i<b+£, {Yi)a<i<b = {yi)a<i<b) 



ie{a,a+l,...,b+l}\{0} 



Furthermore 



IP((Aj)a<j<6+£ = 
-1 



(Xi)a<i<b+e, (Yi) a <i<b — {Vi)a<i<b) 



n 



k=l 



If we compute the sum S in the order: x a , • • • , x-\ first (simplifications using condition ( 19 )) then 
Xb+e,Xb+e-i, ■ ■ ■ ,xe (simplifications using condition p0|)), and finally x\, . . . ,X£-\, we obtain 
eventually: S = h p [xq] Y\ b i=a fi P [yi\- □ 



Corollary 5.9. If both conditions (19) and (20) are satisfied, then every line of the stationary 
space-time diagram consists of i.i.d. random variables. In particular, any two different random 
variables are independent. 

If the neighborhood is M = {0, 1}, the spatial properties of Section [4] remain for a general 
alphabet (existence of transversal PCA, properties of triangles,...). For other neighborhoods, 
there is no natural transversal PCA. 
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6 Cellular automata 



A cellular automaton (CA) is a PCA in which the transition function / is such that, for all 
x G , the probability measure f(x) is concentrated on a single letter of the alphabet. Thus, 
the transition function of a CA can be described by a mapping / : A M — > A, and the CA can 
be viewed as a deterministic mapping F : A 1, — > A z . 

Cellular automata are classical and relevant mathematical objects: they are precisely the map- 
pings from A z to A 1, which are continuous (with respect to the product topology) and commute 
with the shift, see [6]. 



6.1 Known results 

Definition 6.1. A cellular automaton of transition function f : A^ — > A, where the neigh- 
borhood is of the form N = {£,..., r — 1, r} for some i < r, is left-permutative (resp. 
right-permutativej if, for all w = W( • • • w r -i G A r ~ t , the mapping from A to A defined by: 
a i — > f(aw) (resp. a i — > f(wa)), is bijective. A CA is permutative if it is either left or 
right-permutative. 



Let F : A — > A be a permutative CA. The existence of the bijections, see Definition 6.1 has 
two direct consequences: (i) F is surjective; (ii) the uniform measure is invariant: H1/2F = 
In fact, these last two properties are equivalent. 

Proposition 6.2 (Hedlund [B]). Let F be a cellular automaton. We have: 

F is surjective -<=^> ^1/2-^ = ^1/2 • 



There exist surjective CA which are non permutative. Consider, for instance, the mapping 
F : {0, 1} Z — ► {0, 1} Z , defined as follows. Set A = 10010 and B = 11000. Observe that 
the two patterns A and B do not overlap. From a configuration u G {0, 1} Z , we get its image 
Fq(u) by changing each occurence of A into B, resp. of B into A. Clearly, the mapping Fq 
can be defined as a cellular automaton with neighborhood N = {—4, , . . . , 0, . . . , 4}. Also, F is 
surjective but not permutative. 



Let us present a recent result which refines Proposition 6.2 Given a finite and non-empty word 
u G A + , let u 1, = ■■■uuu--- G A z be a periodic bi-infinite word of period u (the starting 
position is indifferent). If F : A z — > A 1, is a CA, then F(u z ) = v z for some word v with 
\v\ = \u\. For simplicity, we write v = F(u). 

Theorem 6.3 (Kari-Taati [7]). Consider a CA F on the alphabet A. The Bernoulli product 
measure fi p , p = (pi)ieA> Pi > for all i, is invariant for F if and only if: 

(i) F is surjective and (ii) Vu G A + , |it|jlog(pi) = \F (u)\i\og(pi) . 

ieA ieA 

Let us mention two consequences of the above results. 

If a cellular automaton has an invariant Bernoulli product measure ji p (pi > for all i), then 
the uniform measure is also invariant. 

A cellular automaton F is number- conserving if: Vu G A + , Vi G A, \u\i = \F(u)\i. A surjective 
and number-conserving CA admits all Bernoulli product measures fi p as invariant measures. 
For instance, the CA Fq, defined above, is surjective and number-conserving. Therefore, all the 
Bernoulli product measures are invariant for Fq. 
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6.2 Link with the conditions for PCA 



The results in Sections 3 4|5 give conditions for a PCA to admit invariant Bernoulli product 
measures. The above results, Section 6.1, give conditions for a CA to admit invariant Bernoulli 
product measures. The natural question is whether we obtain the latter conditions by special- 
izing the former ones. 



Recall that the conditions (jl9h or (20) of Theorem 5.6 are sufficient for the Bernoulli product 



measure \i v (Vi G A,pi > 0) to be invariant for the PCA F. Let us specialize these conditions 
to cellular automata, that is, let us assume that all the coefficients 9„ _ _• are equal to or 1. 

X(J L £— 1 t 



Lemma 6.4. A cellular automaton satisfies condition (19), resp. (20), if and only if it is 
right-permutative, resp. left-permutative. 



Proof. Consider a CA (transition function /) satisfying condition (19) for some p = (pi)i£A- 
Set J = {j E A | Pj = mm ieA pi} and consider j £ J. The equality pj = J2ieAP* ' ^ -^_ii> 
together with the constraints G J Xo ... Xe 1 j £ {Oil}) implies that there must be exactly one index 
k £ J such that O 3 mX k = 1, i.e. f(xo, ■ ■ ■ , a^-i, k) = j. By repeating the argument, we 
obtain that for all xq • • ■ sc£_i, the mapping j h-> J(xq, . . . , xi—x,j) restricted to J is a bijection. 
We now proceed by considering the set of indices J2 = {j G A — J \ Pj = min ie _4\ jPi}, and so 
on. □ 



To summarize, we recover the permutative CA. On the other hand, the sujective but non- 



permutative CA are not captured by the sufficient conditions of Theorem 5.6 



is 



For a left-permutative CA (resp. right-permutative), the transversal CA, see Section 3.3 
right-permutative (resp. left-permutative), and explicitly computable. Moreover, it is well- 
defined even if the space-time diagram is not assumed to be stationary. We recover here a folk 
result. 

In the special case A = {0, 1} and M = {0, 1}, all the surjective CA are permutative. So in 
this case, we recover all the surjective CA. This is consistent with the fact that in this case, the 



conditions of Theorem 5.6 are necessary and sufficient (see Theorem 3.2). 



Remark. Condition (19) can be interpreted as "being right-permutative in expectation" for a 



PCA. And similarly, condition (20) amounts to "being left-permutative in expectation" . 



7 Related open issues 

Consider a PCA of alphabet and neighborhood of size 2. Under the relations Q or ( |10[ ), it 
has an explicit invariant measure with a simple form (Bernoulli product or Markovian). The 
conditions Q and (10) are of codimension 1 in the parameter space. What happens for other 



values of the parameters? Is it still possible to give an explicit description of the invariant 
measure? This is an open question. It has been deeply investigated for the family of PCA's 
defined by: 6>oo = #01 = #10 = a, 0n = 1 — a, for some a S (0, 1). Observe that neither Q nor 



(10) is satisfied except in the trivial case a = 1/2. The specific interest for these PCA is due 
to a connection with directed animals and percolation theory first noticed by Dhar pE], see also 
[21 [8]. More specifically, determining explicitly the invariant measure for the above PCA would 
enable to: (1) compute the area and perimeter generating function of directed animals in the 
square lattice; (2) compute the directed site-percolation threshold in the square lattice. The 
most recent efforts to compute the invariant measure can be found in |9]. 
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